function y=Xap(x)
    y=exp(-3.837)*0.602*x^(0.602-1);
end